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^ ' 1 Introduction 

^ I These notes present a description of quantum chaotic eigenstates, tliat is bound 

W ■ states of a quantum dynamical system, the classical limit of which is chaotic. The 

^ . classical dynamical systems we will be dealing with are mostly of two types: geodesic 

c^ I flows on euclidean domains ("billiards") or compact riemannian manifolds, and 

canonical transformations on a compact phase space; the common feature is the 
"chaoticity" of the dynamics. The corresponding quantum systems will always be 
^ ! considered in the semiclassical (or high-frequency) regime, in order to establish a link 

>• ' with the classical dynamics. As a first illustration, we plot below two eigenstates of a 

paradigmatic system, the Laplacian on the stadium billiard, with Dirichlet boundary 
conditional!. 

</^ I The study of chaotic eigenstates makes up a large part of the field of quan- 

lO ■ tum chaos. It is somewhat complementary with the contribution of J.Keating, who 

O I will focus on the statistical properties of quantum spectra, another major topic in 

quantum chaos. I do not include the study of eigenstates of quantum graphs (a re- 
cent interesting development in the field), since this question should be addressed 
in U.Smilansky's lecture. Although these notes are purely theoretical, H.-J. Stock- 
r> I mann's lecture will show that the questions raised here have direct experimental 

c^ ■ applications (his lecture should present experimentally measured eigenmodes of 2- 

and 3-dimensional "quantum billiards"). 

One common feature of the chaotic eigenf unctions (except in some very specific 
systems) is the absence of explicit (or even approximate) formulas for them. One 
then has to resort to indirect, rather unprecise approaches to describe these eigen- 
states. We will use various analytic tools or points of view: deterministic/statistical, 
macro/microscopic, pointwise/global properties, generic/specific systems. The level 
of rigour in the results varies from mathematical proofs to heuristics, generally sup- 
ported by numerical experiments. The necessary selection of results reflects my per- 
sonal view or knowledge of the subject, it omits several important developments, 
and is more "historical" than sharply up-to-date. The list of references is qisw but 
in no way exhaustive. 
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*I am grateful to E.Bogomolny, who allowed me to reproduce several plots from 1261 . The author has been 
partially supported by the Agence Nationale de la Recherche under the grant ANR-09-JCJC-0099-01. These notes 
were written while he was visiting the Institute of Advanced Study in Princeton, supported by the National Science 
Foundation under agreement No. DMS-0635607. 

^The eigenfunctions of the stadium plotted in this article were computed using a code nicely provided to me by 
E. Vergini, which uses the scaling method invented in |100| . 
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Figure 1: Two eigenfunctions of the Dirichlet Laplacian in the stadium bilhard, with wavevectors 
k — 60.196 and k = 60.220 (see ([9])). Large values of |'!/'(a:)p correspond to dark regions, while 
nodal lines are white. While the left eigenfunction looks relatively "equidistributed" , the right one 
is scarred by two symmetric periodic orbits (see 



These notes are organized as follows. We introduce in section |2] the classical 
dynamical systems we will focus on (mostly geodesic flows and maps on the 2- 
dimensional torus), mentioning their degree of "chaos". We also sketch the quan- 
tization procedures leading to quantum Hamiltonians or propagators, the eigen- 
states of which we want to understand. We also mention some properties of the 
semiclassical/high-frequency limit. In section |3] we describe the macroscopic prop- 
erties of the eigenstates in the semiclassical limit, embodied by their semiclassical 
measures. These properties include the quantum ergodicity property, which for some 
systems (with arithmetic symmetries) can be improved to quantum unique ergodic- 
ity, namely the fact that all high-frequency eigenstates are "flat" at the macroscopic 
scale; on the opposite, some specific systems allow the presence of exceptionally local- 
ized eigenstates. In section H] we focus on more refined properties of the eigenstates, 
many of statistical nature (value distribution, correlation functions). Very little is 
known rigorously, so one has to resort to heuristical models of random wavefunc- 
tions to describe these statistical properties. The large values of the wavefunctions or 
Husimi densities are discussed, including the scar phenomenon. Section \5\ discusses 
the most "quantum" or microscopic aspect of the eigenstates, namely their nodal 
sets, both in position and phase space (Husimi) representations. Here as well, the 
random state models are helpful, and lead to interesting questions in probability 
theory. 

2 What is a quantum chaotic eigenstate? 

In this section we first present a general definition of the notion of "chaotic eigen- 
state". We then focus our attention to geodesic flows in euclidean domains or on 
compact riemannian manifolds of negative curvature, which form the simplest sys- 
tems proved to be chaotic. Finally we present some discrete time dynamics (chaotic 
canonical maps on the 2-dimensional torus). 



2.1 A short review of quantum mechanics 

Let us start by recalling that classical mechanics on the phase space T*R°' can 
be defined, in the Hamiltonian formalism, by a real valued function H{x,p) on 
that phase space, called the Hamiltonian. We will always assume the system to be 
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autonomous, namely the function H to be independent of time. This function then 
generates the flo\\o 

x(t) =^(x(t),p(t)) = <l>*j^(x(0)), tGM, x(t)GT*M^ 

by solving Hamilton's equations (as usual / is the time derivative): 

ciM flFf 

i,W = ^(x(t)), Pj{t) = - — {^{t)), j = l,...,d. (1) 

This flow preserves the symplectic form 'Yl,jdpj A dxj, and the energy shells £e = 
H~\E). 

The corresponding quantum mechanical system is defined by an operator Hji 
acting on the (quantum) Hilbert space l-L = L'^{M.'^,dx). This operator can be for- 
mally obtained by replacing coordinates x,p by operators: 

Hn = Hixn,Pfi), (2) 

where xn is the operator of multiplication by x, while the momentum operator 
Pfi = ^V is conjugate to x through the h-Fouher transform J^fj. The notation ([2]) 
assumes that one has selected a certain ordering between the operators xn. and pri, in 
physics one usually chooses the fully symmetric ordering, also called the Weyl quan- 
tization: it has the advantage to make Hfi formally self-adjoint on L^(M'^). Quan- 
tization procedures can also be defined when the euclidean space M*^ is replaced 
by a compact manifold M, by gluing together quantizations defined on different 
coordinate patches. 

The quantum dynamics, which governs the evolution of the wavefunction ilj{t) G 
"H describing the system, is then given by the Schrodinger equation: 

^^^=[Hnm). (3) 

Solving this linear equation produces the propagator, that is the family of unitary 
operators on L'^{M.'^), 

Ul = exp{-iHn/h), teR. 



Remark 1 In physical systems, Planck's constant h is a fixed number, which is of 
order 10"'^^ in standard units. However, if the system (atom, molecule, "quantum 
dot") is itself microscopic, the value of h may be comparable with the typical action 
of the system, in which case it is more natural to select units in which h = 1. Our 
point of view throughout this work will be the opposite: we will assume that h is 
(very) small compared with the typical action of the system, and most results will be 
valid asymptotically, in the semiclassical limit /i — )■ 0. 

2.2 Quantum-classical correspondence 

At this point, let us introduce the crucial semiclassical property of the quantum evo- 
lution: it is called (in the physics literature) the quantum-classical correspondence. 



We always assume that the flow is complete, that is it does not blow up in finite time. 
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while in mathematics this result is known as Egorov's theorem. This property states 
that the evolution of observables approximately commutes with quantization. For 
us, a (classical) observable is a smooth, compactly supported function on phase 
space / e C^(T*]R'^), while the corresponding quantum observable fn, is obtained 
through the above mentioned quantization procedure. The evolution of classical and 
quantum evolutions are defined by duality with that of particles/wavefunctions: 

/(t) = /o$^, f,{t) = U^'fnUl 

The quantum-classical correspondence connects these two evolutions: 

Vt e M, hit) = fit)n + OL2^L<e^^'^h) , (4) 

where the exponent F > depends on the instability of the classical flow. 

The most common form of dynamics is the motion of a scalar particle in an 
electric potential V{x). It corresponds to the Hamiltonian 

H{x,p) = h qV{x), quantized into Hn = \- qV{x) . (5) 

2m 2m 

We will usually scale the mass and electric charge to m, = q = 1, keeping h small. 
Since the Hamilton flow ([1]) leaves each energy energy shell Se invariant, we may 
restrict our attention to the flow on a single shell. We will be interested in cases 
where 

1. the energy shell Se is bounded in phase space (that is, both the positions and 
momenta of the particles remain finite at all times). This is the case if V{x) is 
confining (y{x) — t- oo as |x| — )■ oo). 

2. the flow on Se is chaotic (and this is also the case on the neighbouring shells 
Se+e)- "Chaos" is a vague word, which we will make more precise below. 

The first condition implies that, provided h is small enough, the spectrum of H/i is 
purely discrete near the energy E, with eigenstates ipnj £ L'^iM.'^) {bound states). 
Besides, fixing some small e > and letting ^ — )■ 0, the number of eigenvalues 
Efij E [E — e, E + e] typically grows like Ch~'^. Under the second condition, the 
eigenstates with energies in this interval can be called quantum chaotic eigenstates. 
Below we describe several degrees of "chaos" , which regard the long time prop- 
erties of the classical flow. These properties are relevant when describing the eigen- 
states of the quantum system, which form the "backbone" of the long time quantum 
dynamics. The main objective of quantum chaos consists in connecting, in a precise 
way, the classical and quantum long time (or time independent) properties. 

2.3 Various levels of chaos 

For most Hamiltonians of the form ([5]) (e.g. the physically relevant case of a hydro- 
gen atom in a constant magnetic field), the classical dynamics on bounded energy 
shells Se involves both regular and chaotic regions of phase space; one then speaks 
of a mixed dynamics on Se- The regular region is composed of a number of "islands 
of stability" , made of quasiperiodic motion structured around stable periodic orbits; 
these islands are embedded in a "chaotic sea" where trajectories are unstable (they 
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have a positive Lyapunov exponent). These notions of "island of stabihty" versus 
"chaotic sea" are rather poorly understood mathematically, but have received com- 
pelling numerical evidence [TTj. The main conjecture concerning the corresponding 
quantum system, is that most eigenstates are either localized in the regular region, 
or in the chaotic sea [83]. To my knowledge this conjecture remains fully open at 
present, in part due to our lack of understanding of the classical dynamics. 

For this reason, I will restrict myself (as most researches in quantum chaos do) 
to the case of systems admitting a purely chaotic dynamics ouSe- I will allow various 
degrees of chaos, the minimal assumption being the ergodicity of the flow $^ on Se, 
with respect to the natural (Liouville) measure on Se- This assumption means that, 
for almost every initial point :s.q ^ £e, the time averages of any observable / converge 
to its phase space average: 

r 1 r ,r^t r ^^. f r, ,r . .def //(X)(5(i/(X)-E)dx 

T^oo 2T J_rp J^^ J 6{H{x} - E) dx 

A stronger chaotic property is the mixing property, or decay of time correlations 
between two observables f,g: 

CfAt) = I 9y<{fo<l>')dfiL- I fdfiLJ gdfiL^^^O. (7) 

The decay rate of Cf^g{t) (or mixing rate) depends on both the flow $* and the 
regularity of the observables /, g. For very chaotic flows (Anosov flows, see §2.4.2p 
and smooth observables, the decay is exponential. 

2.4 Geometric quantum chaos 

In this section we give explicit examples of chaotic flows, namely the geodesic flows 
in certain euclidean billiards compact riemannian manifolds. The dynamics is then 
induced by the geometry, rather than by a potential. Both the classical and quantum 
properties of these systems have been investigated a lot in the past 30 years. 

2.4.1 Billiards 

The simplest form of ergodic system occurs when the potential V{x) is an infinite 
barrier delimiting a bounded domain fi C M'' (say, with piecewise smooth boundary), 
so that the particle moves freely inside Q and bounces specularly at the boundaries. 
For obvious reasons, such a system is called a euclidean billiard. All positive energy 
shells are equivalent to one another, up to a rescaling of the velocity, so we may 
restrict our attention to the shell £ = {{x,p), x E fl, \p\ = 1}. The long time 
dynamical properties only depend on the shape of the domain. For instance, in 2 
dimensions, rectangular, circular or elliptic billiards lead to integrable dynamics: 
the flow admits two independent integrals of motion — in the case of the circle, 
the energy and the angular momentum. A convex billiard with a smooth boundary 
will always admit "whispering gallery" stable orbits, so such billiards cannot be 
fully chaotic. On the opposite, the famous stadium billiard (see Fig. [2]) was proved 
to be ergodic by Bunimovich [31]. Historically, the first euclidean billiard proved 
to be ergodic was the Sinai billiard, composed of one or several circular obstacles 
inside a square (or torus) [97]. These billiards also have positive Lyapunov exponents 
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(meaning that almost all trajectories are exponentially unstable, see the left part of 
Fig. |2]). It has been shown more recently that these billiards are mixing, but with 
correlations decaying at polynomial or subexponential rates [361 EH CS] • 




Figure 2: Left: two trajectories in the stadium billiard, initially very close to one another, and 
then diverging fast (the red arrow shows the initial point). Right: long evolution of one of the 
trajectories. 

The quantization of the broken geodesic flow inside Q is the (semiclassical) 
Laplacian with Dirichlet boundary conditions: 



Hn 



h'A, 



(8) 



Obviously, the parameter h only amounts to a rescaling of the spectrum: an eigen- 
state ipn of ([8]) with energy Efi ^ 1/2 is also an eigenstate of — Af^ with eigenvalue 
k"^ ^ h~'^. Hence, h represents the wavelength of il^n, the inverse of its wavevector k. 
Fixing E = 1/2 and taking the semiclassical limit /i — )■ is equivalent with studying 
the high-frequency or high- wavevector spectrum of — A^. 

The system (|8]) is often called a quantum billiard, although this operator is 
not only relevant in quantum mechanics, but in all sorts of wave mechanics (see 
H.-J. Stockmann's lecture). Indeed, the scalar Helmholtz equation 



A^. + kU, = 



(9) 



may describe stationary acoustic waves in a cavity, and is also relevant to describe 
electromagnetic waves in a quasi-2D cavity, provided one is allowed to separate the 
different polarization components of the electric field. 

Euclidean billiards thus form the simplest realistic quantized chaotic systems, for 
which the classical dynamics is well understood at the mathematical level. Besides, 
the spectrum of the Dirichlet Laplacian can be numerically computed up to large 
values of k using methods specific to the euclidean geometry, like the scaling method 
|100j . For these reasons, these billiards have become a paradigm of quantum chaos 
studies. 



2.4.2 Anosov geodesic flows 

The strongest form of chaos occus in systems (maps or flows) with the Anosov 
property, also called uniformly hyperbolic systems [5]. The first (and main) example 
of an Anosov flow is given by the geodesic flow on a compact riemannian manifold 
{M,g) of negative curvature, generated by the free particle Hamiltonian H{x,p) = 
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|p|^/2. Uniform hyperbolicity — which is induced by the negative curvature of the 
manifold — means that at each point x G £^ the tangent space Ty^£ sphts into the 
vector Xx generating the flow, the unstable subspace E^ and the stable subspace 
E~. The stable (resp. unstable) subspace is defined by the property that the flow 
contracts vectors exponentially in the future (resp. in the past), see Fig. O for some 
C > 0, A > independent of a:; G £^, 



Vt; e E^, Vt > 0, 



\\d^fv\\ <Ce-^'\\v\\. 



(10) 



Anosov systems present the strongest form of chaos, but their ergodic properties 
are (paradoxically) better understood than for the billiards of the previous section. 
The flow has a positive complexity, reflected in the exponential proliferation of long 
periodic geodesies. For this reason, this geometric model has been at the center of 
the mathematical investigations of quantum chaos, in spite of its minor physical 
relevance. Generalizing the case of the euclidean billiards, the quantization of the 





Figure 3: Top left: geodesic flow on a surface of negative curvature (such a surface has a genus 
> 2). Right: fundamental domain for an "octagon" surface r\IHI of constant negative curvature (the 
figure is due to C. McMullen) Botton left: a phase space trajectory and two nearby trajectories 
approaching it in the future or past. The stable/unstable directions at x and $*(x) are shown. The 
red lines feature the expansion along the unstable direction, measured by the unstable Jacobian 

jr(x) = det(d$* r^+)- 



geodesic flow on (M, g) is given by the (semiclassical) Laplace- Beltrami operator 



Hn 



h^X 



'11^ 



acting on the Hilbert space L'^{M,dg) associated with the Lebesgue measure. The 
eigenstates of Hfi with eigenvalues ~ 1/2 (equivalently, the high-frequency eigen- 
states of —Ag) constitute a class of quantum chaotic eigenstates, the study of which 
is not impeded by boundary problems present in billiards. 

The spectral properties of this Laplacian have interested mathematicians work- 
ing in riemannian geometry, PDEs, analytic number theory, representation theory, 
for at least a century, while the specific "quantum chaotic" aspects have emerged 
only in the last 30 years. 
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The first example of a manifold with negative curvature is the Poincare half- 
space EI with its hyperbolic metric ^ ^ ^ , on which the group 5X2 (M) acts iso- 
metrically by Moebius transformations (H is often represented by the equivalent 
hyperbolic disk, like in Fig. |3]). For certain discrete subgroups F of PSL2(M.) (called 
co-compact lattices), the quotient M = F\EI is a smooth compact surface. This 
group structure provides detailed information on the spectrum of the Laplacian (for 
instance, the Selberg trace formula exphcitly connects the spectrum of A^ with the 
periodic geodesies on M, which are themselves represented by the conjugacy classes 
ofF). 

Furthermore, for some of these discrete subgroups F (called arithmetic) , one can 
construct a commutative algebra of Hecke operators on L'^{M), which also commute 
with the Laplacian; it then make sense to study in priority the joint eigenstates of A 
and of these Hecke operators, which we will call the Hecke eigenstates. This arith- 
metic structure allows to obtain extra nontrivial information on these eigenstates 
(see §3.3.3p . so the latter will appear several times along these notes. Their study 
composes a part of arithmetic quantum chaos, a lively field of research. 



2.5 Classical and quantum chaotic maps 

Beside the Hamiltonian or geodesic flows, another model system has attracted much 
attention in the dynamical systems community: chaotic maps on some compact 
phase space V. Instead of a flow, the dynamics is given by a discrete time transfor- 
mation K, : V -^ V. Because we want to quantize these maps, we require the phase 
space V to have a symplectic structure, and the map k to preserve this structure (in 
other words, k is an invertible canonical transformation on V). The advantage of 




p=sm((p) 





Figure 4: Top: Poincare section and the associated return map constructed from a Hamiltonian 
flow on £e- Bottom: boundary map associated with the stadium bilhard. 



studying maps instead of flows is multifold. Firstly, a map can be easily constructed 
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from a flow by considering a Poincare section S transversal to the flow; the induced 
return map k^ : T, ^ Tj, together with the return time, contain all the dynamical in- 
formation on the flow. Ergodic properties of chaotic maps are usually easier to study 
than their flow counterpart. For billiards, the natural Poincare map to consider is 
the boundary map ks defined on the phase space associated with the boundary, 
T*dQ. The ergodic properties of this boundary map were understood, and used to 
address the case of the billiard flow itself |14j . 

Secondly, simple chaotic maps can be defined on low- dimensional phase spaces, 
the most famous ones being the hyperbolic symplectomorphisms on the 2-dimensional 

torus. These are defined by the action of a matrix "S* = j , I with integer entries, 

determinant unity and trace a + d > 2 (equivalently, 5* is unimodular and hyper- 
bolic). Such a matrix obviously acts on x = {x,p) G T^ linearly, through 



i^s{^) 



[ax 



bp, ex + dp) mod 1 . 



(12) 



A =chemaac v,ew of .. fo. the famous ArnoUS cat ».p S., = (j ^ . displayed 

in Fig. [5l The hyperbolicity condition implies that the eigenvalues of S are of the 
form {e^'^} for some A > 0. As a result, ks has the Anosov property: at each point x, 
the tangent space TxT^ splits into stable and unstable subspaces, identified with the 
eigenspaces of 5, and ±A are the Lyapunov exponents. Many dynamical properties of 
Kg can be explicitly computed. For instance, every rational point x G T^ is periodic, 
and the number of periodic orbits of period < n grows like e'^" (thus A also measures 
the complexity of the map). This linearity also results in the fact that the decay 
of correlations (for smooth observables) is superexponential. To obtain a generie 








P' 


I 




.-/"V' 






V- 


/ 


"'■;/ 


^ ^^ 




,'' 




^''' 


'"/'' "" 






^ ■' 


V"' 


,-'' 


X 


■^ 


^-' 


^ 




\ 




/''' 








\ 





OnT 




Figure 5: Construction of Arnold's cat map ks^^j on the 2-torus, obtained by periodizing the Unear 
transformation on R^. The stable/unstable directions are shown (kindly provided by F. Faure). 

Anosov diffeomorphism of the 2-torus, one can smoothly perturb the linear map Ks'- 
given a Hamiltonian H G C°°(T^), the composed map $^ o ks remains Anosov if 
e is small enough, due to the structural stability of Anosov diffeomorphisms. For 
any such Anosov diffeomorphism, the decay of correlations is exponential. Another 
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Figure 6: Schematic view of the baker's map ([T3|) . The arrows show the contraction/expansion 
directions. 

family of canonical maps on the torus was also much investigated, namely the so- 
called baker's maps, which are piecewise linear. The simplest (symmetric) baker's 
map is defined by 



kb{x,p) 



(2xmodl,|), 
{2x mod 1, 2±i; 



< X < 1/2, 
1/2 < X < 1. 



(13) 



This map is conjugate to a very simple symbolic dynam,ics, namely the shift on 
two symbols. Consider the binary expansions of the coordinates x = 0.aia2 ■ ■ ■ , 
p = 0.(3i/32 ■ ■ ■ . Then the map {x,p) i— )■ kb{x,p) is equivalent with the shift to the 
left on the bi- infinite sequence ■ ■ ■ /32/3i ■ «ia2 ■ ■ ■ • This conjugacy allows to easily 
identify all periodic orbits, construct a large set of nontrivial invariant probability 
measures, and prove that the map is ergodic and mixing w.r.t. these measures. All 
trajectories not meeting the discontinuity lines are uniformly hyperbolic. 

Simple canonical maps have also be defined on the 2-sphere phase space (like the 
kicked top), but their chaotic properties have, to my knowledge, not been rigorously 
proven. Their quantization has been intensively investigated, mostly numerically 



2.5.1 Quantum maps on the 2-dimensional torus 

As opposed to the case of Hamiltonian flows, there is no natural rule to quan- 
tize a canonical map on a compact phase space V. Already, associating a quantum 
Hilbert space to this phase space is not obvious. Therefore, from the very beginning, 
quantum maps have been defined through somewhat arbitrary (or rather, ad hoc) 
procedures, often specific to the considered map k : V ^ V. Still these recipes are 
always required to satisfy a certain number of properties: 

• one needs a sequence of Hilbert spaces {'Hiy)iy^f4 of dimensions N. Here N is 
interpreted as the inverse of Planck's constant, in agreement with the heuristics 
that each quantum, state occupies a volume h!^ in phase space. We also want to 
quantize observables / G C{V) into hermitian operators /at on T-Ln. 

• For each A^ > 1, the quantization of n is given by a unitary propagator Un{k) 
acting on T-Ln- The whole family {UN{'i))N>i is called the quantum map asso- 
ciated with K. 
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• in the semiclassical limit N ^ h^^ ^ oo, this propagator satisfies some form of 
quantum-classical correspondence. Namely, for some (large enough) family of 
observables / on V, we should have 

VnGZ, U^^fNU^ = {f7^^)^ + On{N-') asiV->oo. (14) 

The condition ( !T4|) is the analogue of the Egorov property (j4]) satisfied by the propa- 
gator Uri associated with a quantum Hamiltonian, which quantizes the stroboscopic 
map X H- <|)^(x). 

Let us briefly summarize the explicit construction of the quantizations Un{h), 
for the maps k : T^ — ;■ T^ presented in the previous section. Let us start by con- 
structing the quantum Hilbert space. One can see T^ as the quotient of the phase 
space T*M. = R^ by the discrete tranlations x i— )■ x+n, n G Z,^. Hence, it is natural 
to construct quantum states on T^ by starting from states ip G i5'(]R), and requiring 
the following periodicity properties 

It turns out that these two conditions can be satisfied only ii h = (2ttN)^^, A^ G N, 
and the corresponding distributions then form a vector space Hn of dimension A^. 
A basis of this space is given by the Dirac combs 

^^(^) = ;^ E ^(^ - ^ - '^)' i = 0,...,N-l. (15) 



It is natural to equip T-Ln with the hermitian structure for which the basis {eg, i = 
0, . . . , A^ — 1} is orthonormal. The components of a vector ip G "Hat in this "position 
basis" will be denoted by ip^i/N) = {ee,ip). 

Let us now explain how the symplectomorpisms k^ are quantized [51]. Given 
a unimodular matrix S, its action on M^ can be generated by a quadratic polyno- 
mial Hs{x,p); this action can thus be quantized into the unitary operator Uri{S) = 
exp{—iHsfi/fi) on L^(R). This operator also acts on distributions iS'(M), and in par- 
ticular on the finite subspace "Hat. Provided the matrix S satisfies some "checker- 
board condition", one can show (using group theory) that the action of Uji{S) on 
"Htv preserves that space, and acts on it through a unitary matrix Un{S). The fam- 
ily of matrices {Un{S))n>i defines the quantization of the map ks on T^. Group 
theory also implies that an exact quantum-classical correspondence holds (that is, 
the remainder term in flT^ vanishes) , and has other important consequences regard- 
ing the operators Un = Un{i^s) (for each A^ the matrix Un is periodic, of period 
Tn < 2A^). Explicit expressions for the matrices Un{S) can be worked out, they 
depends sensitively on the factorization into primes of the integer A^. 

The construction of the quantized baker's map (1131) proceeds differently. An 
Ansatz was proposed by Balasz-Voros [13], which takes the following form in the 
basis (HgH: 

UnM = F*^ (^^/^ ^^^^) , (16) 

where F/v is the A^- dimensional discrete Fourier transform. This Ansatz is obviously 
unitary. It was obtained by discretizing the phases e"^^^''^'-*/'*, where the function 

^We assume that A'^ is an even integer 
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(p{p',x) locally generatecl the map ( 1T3|) . A proof that the matrices Un{i^b) satisfy 
the Egorov property (iMl) was given in [4lj. 

Once we have constructed the matrices Un{i^) associated with a chaotic map k, 
their eigenstates {ipNj, j = 1, • • • , N} enjoy the role of quantum chaotic eigenstates. 
They are of quite different nature from the eigenstates of the Laplacian on a manifold 
or a billiard: while the latter belong to L'^{M) or L^(f2) (and are actually smooth 
functions), the eigenstates ipN,j are iV-dimensional vectors. Still, part of "quantum 
chaos" has consisted in developing common tools to analyze these eigenstates, in 
spite of the different functional settings. 

3 Macroscopic description of the eigenstates 

In this section we study the macroscopic localization properties of chaotic eigen- 
states. Most of the results are mathematically rigorous. In the case of the semiclas- 
sical Laplacian ([8]) on a billiard Q we ask the following question: 

Consider a sequence {il)ji)h^o of normalized eigenstates of Hn^ with energies 
Ef^ ^ 1/2. For A (ZVL a, fixed subdomain, what is the probability that the 
particle described by the stationary state ipji lies inside A! How do the 
probability weights /^ \ipfi{x)\'^ dx behave when ^ — )■ 0? 

This question is quite natural, when contemplating eigenstate plots like in Fig. [H 
Here, by macroscopic we mean that the domain A is kept fixed while /i — )■ 0. 

One can obviously generalize the question to integrals of the type /^ /(x) |^;i(a;) p dx, 
with /(x) a continuous test function on Vt. Assuming the volume of Vt is normal- 
ized, this integral can be interpreted as the quantum average^ {ijjfi, fn^'h), where the 
quantum observable ffi is just the multiplication operator by /(x). It proves useful 
to extend the question to phase space observables /(x,pjj: what is the behaviour of 
the quantum averages 

<(/) = (V^/.,A^ft), feC^{T*n), in the limit n^O? (17) 

Since the quantization procedure / i— )■ /;» is linear, these averages define a distribu- 
tion fi^^ = ly^^(x) (ix on T*Q, called the Wigner distribution of the state ipfi (the 
density IV^^(x) is called the Wigner function). Although this distribution is gener- 
ally not positive, it is interpreted as a quasi-probability density describing the state 
iph in phase space. 

On the euclidean space, one can define a nonnegative phase space density asso- 
ciated to the state iIjji: the Husimi measure (and function): 



/xj^ = n^.ip.) d^, n^.i^) = (27rn)-'^/2 |(y,,,^,)|2 



where the Gaussian wavepackets (y^x G L'^iM.'^) are defined in (pT]) . This measure can 
also be obtained by convolution of jj,^ with the Gaussian kernel e~'^~^' /''. In these 
notes our phase space plots show Husimi functions. 



*This means that, in some region, {x',p') = Kg(x,p) is defined by solving x' = d„/ip{p',x), p = dx<fi{p' ,x). 

^These are also referred to as the (diagonal) matrix elements for the "matrix" /fi. 

^If n C M'* is replaced by a compact riemannian manifold M, we use some quantization scheme f ^^ fa- 
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These questions lead us to the notion of phase space locahzation, or microlo- 
calizatioixX We will say that the family of states {ijjn) is microlocalized inside a set 
B C T*VL if, for any smooth observable f{x,p) vanishing near B, the quantum 

averages {ipn, fn'iph) decrease faster than any power of h when ^ — )■ 0. 

Microlocal properties are not easy to guess from plots of the spatial density 
\ipj (x) p like in Fig. [H but they are natural to study if we want to connect quantum to 
classical mechanics, since the latter takes place in phase space rather than in position 
space. Indeed, the major tool we will use is the quantum-classical correspondence 
(jl]); for all the flows we consider, a generic spatial test function f{x) evolves through 
$* into a genuine phase space function ft{x,p). Microlocal properties are easier to 
visualize on 2-dimensional phase spaces (see below the figures on the 2-torus). 

3.1 The case of completely integrable systems 

In order to motivate our further discussion of chaotic eigenstates, let us first recall 
a few facts about the antipodal systems, namely completely integrable Hamiltonian 
flows. For such systems, the energy shell Se is foliated by d- dimensional invariant 
Lagrangian tori. Each such torus is characterized by the values of d independent 
invariant actions /i, . . . , 7^, so let us call such a torus Tf. The WKB theory allows one 
to explicitly construct, in the semiclassical limit, precise quasimodes of Hri associated 
with some of these toru, that is normalized states il)f = ipj-^ f satisfying 

Hn^Pf=Efi;j-+0{h'^), (19) 

with energies Ef ^ E. Such a Lagrangian (or WKB) state ipf takes the following 
form (away from caustics): 

L 

^/(^) = H M^-. K) eM^Si{x)/h) . (20) 

(.=1 

Here the functions Si{x) are (local) generating functional for Tf, and each Ai{x] K) = 
y4°(a;) + hA\{x) + ■ ■ ■ is a smooth amplitude. 

From this explicit expression one can check that the state ipf is microlocalized on 
Tf. On the other hand, our knowledge of t/^^i is much more precise than the latter fact; 
indeed, one can easily construct states microlocalized on Tf, which are very different 
from the Lagrangian states ipf. For instance, for any point xq = (xo,Po) ^ Tf the 
Gaussian wavepacket (or coherent state) 

<^xo(a;) = (7rfi)-i/^e-l"-"ol'/2?^e*Po-"/^' (21) 

is microlocalized on the single point xq, and therefore on Tf. This example just 
reflects the fact that a statement about microlocalization of a sequence of states 
provides much less information than a formula like (120|) . 



^The prefix micro mustn't mislead us: we are still dealing with macroscopic localization properties of ?/)ft! 
^The "quantizable" tori T? satisfy Bohr-Sommerfeld conditions li = 2-iTh(ni + ctj), with Ui a "L arbitrary and 
cii a [0, 1] fixed Maslov indices. 

^Above some neighbourhood U G W^, the torus T? is the union of L lagrangian leaves {{x,'VSi(x)), x S U}, 

e=i,...,L 
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Figure 7: Left: one eigenmode of the circle billiard, microlocalized on a torus Tj- (the invariants 
/i,/2 can be taken to be the energy and the angular momentum). Center: Husimi density of a 
simple Lagrangian state on T^, a momentum eigenstate microlocalized on the lagrangian {■C = ■fo}- 
Right: Husimi density of the Gaussian wavepacket (HI]). Stars denote the zeros of the Husimi 
density, (reprinted from [82]) 



3.2 Quantum ergodicity 

In tlie case of a fully chaotic system, there is no explicit formula describing the eigen- 
states, or even the quasimodes, of Hfi. However, macroscopic informations on the 
eigenstates can be obtained indirectly, using the quantum-classical correspondence. 
The main result on this question is a quantum analogue of the ergodicity property 
(P) of the classical flow, and it is a consequence of this property. For this reason, it 
has been named Quantum Ergodicity by Zelditch. Loosely speaking, this property 
states that almost all eigenstates ipfi with Ef^ ^ E will become equidistributed on 
the energy shell Se, in the semiclassical limit, provided the classical flow on Se is 
ergodic. We give below the version of the theorem in the case of the Laplacian on a 
compact riemannian manifold, using the notations of (^. 

Theorem 1 [Quantum ergodicity] Assum,e the geodesic flow on (M, g) is ergodic 
w.r.to the Liouville measure. Then, for any orthonormal eigenbasis {'4'j)j>o of the 
Laplacian, there exists a subsequence S C N of density 1 (that is, limj_j.oo j ' = 
1), such that 



v/ E cr(M), 






f{x,p)dfXL{x,p) 



where hl is the normalized Liouville measure on £ = E\i2, and hj = k, ^ the inverse 
wavevector ofipj. 

The statement of this theorem was first given by Schnirelman (using test functions 
fix)) [93], the complete proof was obtained by Zelditch in the case of manifolds 
of constant negative curvature T\M |104j . and the general case was then proved by 
Colin de Verdiere EH. This theorem is "robust": it has been extended to 



quantum ergodic billiards ^\ 1109] 

quantum Hamiltonians H/i, such that the flow $^ is ergodic on Se in some 
energy interval [53] 
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• quantized ergodic diffeomorphisins on the torus [28] or on more general compact 
phase spaces |106j 

• a general framework of C* dynamical systems |105j 

• a family of quantized ergodic maps with discontinuities [75], in particular the 
baker's map [41] 



• certain quantum graphs [T7] 

Let us sketch the proof of the above theorem. We want to study the statistical 
distribution of the matrix elements fifif) ~ ii^jy ffiji^j) i^ ^^^ range {kj < K}, 
with K ^ 1. The first step is to estimate the average of this distribution. It is 
estimated by the generalized Weyl law: 

Eft(/)~^^^|^A"'//*- '<^^- (22) 

kj<K ^ '^' -^^ 

where ad is the volume of the unit ball in W^ [56]. In particular, this asymptotics 
allows to count the number of eigenvalues kj < K: 

#{%</0~^^^/f^ /f^oo. (23) 

and shows that the average of the distribution {/ij^(/), kj < K} converges to the 
phase space average /iL(/) when K ^ oo. 

Now, we want to show that the distribution is concentrated around its average. 
This can be done by estimating its variance 

which is often called the quantum variance. Because the ipj are eigenstates of the 
propagator [/^, we may replace f^. by its time average up to some large time T, 

hfij = ^ / ^h-fhj^hj dt , 

without modifying the matrix elements. Then we apply the Cauchy-Schwarz inequal- 
ity: 

\{ipj, Ai/jj\'^ < {ipj,A* Alp), for any bounded operator A, 

to get the following upper bound for the variance: 

"rri j — i f,.^j^ 

The Egorov theorem (jl]) shows that the product operator on the right hand side is 
approximatly equal to the quantization of the function |/r — /iL(/)p, where fr is 
the time average of the classical observable /. Applying the generalized Weyl law 
(1221) to this function, we get 

Var^(/) < /ii(|/r - /iL(/)n + Ot{K-') , K ^ oo. 
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Finally, the ergodicity of the classical flow imphes that f^iilfT — AiL(/)p) converges 
to zero when T — )■ oo. By taking T, and then K, large enough, the above right hand 
side can be made arbitrary small. This proves that the quantum variance Yaxxif) 
converges to zero when i^ — )■ oo. A standard Chebychev argument is then used 
to extract a dense subsequence such that {ipn, fni'h) — ^ f^L^f), and a flnal diagonal 
extraction procedure yields a dense subsequence enjoying the same property for all 
smooth observables. D 

A more detailed discussion on the quantum variance and the distribution of 
matrix elements {/ij(/), kj < K} will be given in §4.1.3[ 

3.3 Beyond QE: Quantum Unique Ergodicity vs. strong scarring 

3.3.1 Semiclassical measures 

Quantum ergodicity can be conveniently expressed by using the concept of semi- 
classical measure. Remember that, using a duality argument, we associate to each 
eigenstates ipj a Wigner distribution fiY on phase space. The quantum ergodicity 
theorem [1] can be rephrased as follows: 

There exists a density- 1 subsequence 5 C N, such that the sequences of 
Wigner distributions {fiY)j<^s weak-* (or vaguely) converges to the Liou- 
ville measure on S. 

For any compact riemannian manifold {M,g), the sequence of Wigner distributions 
{fiY)jeN remains in a compact set in the weak-* topology, so it is always possible 
to extract an inflnite subsequence (/ij)jG5 vaguely converging to a limit distribution 
fisc, that is 

V/ G C^ {T*M) , lim / / rf/if = f f rf/i,, . 

Such a limit distribution is necessarily a probability measure on £., and is called 
a semiclassical measure of the manifold M. The quantum-classical correspondence 
implies that /igc is mwanant through the geodesic flow: (<l>*)*/isc = ^sc- The semiclas- 
sical measure /i^c represents the asymptotic (macroscopic) phase space distribution 
of the eigenstates {'il)j)j^s- It is the major tool used in the mathematical literature 
on chaotic eigenstates (see below). The deflnition can be obviously generalized to 
any quantized Hamiltonian flow or canonical map. 

3.3.2 Quantum Unique Ergodicity conjecture 

The quantum ergodicity theorem provides an incomplete information, which leads 
to the following question: 

Question 1 Do all eigenstates become equidistrihuted in the semiclassical limit? Equiv- 
alently, is the Liouville measure the unique semiclassical measure for the manifold 
M? On the opposite, are there exceptional subsequences converging to a semiclas- 
sical measure fisc ¥" f^L? 

This question makes sense if the geodesic flow admits invariant measures different 
from fiL (that is, the flow is not uniquely ergodic). Our central example, manifolds 
of negative curvature, admit many different invariant measures, e.g. the singular 
measures //^ supported on each of the (countably many) periodic geodesies 7. 
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This question was already raised in [37], where the author conjectured that 
no subsequence of eigenstates can concentrate along a single periodic geodesic, in 
other words /i^ cannot be a semiclassical measure. Such an unlikely subsequence 
was later called a strong scar by Rudnick and Sarnak ^U\ , in reference to the scars 
discovered by Heller on the stadium billiard (see §4.21) . In the same paper the authors 
formulated a stronger conjecture: 

Conjecture 1 [Quantum unique ergodicity] Let (M, g) be a compact riemannian 
manifold with negative sectional curvature. Then all high-frequency eigenstates of the 
Laplacian become equidistributed with respect to the Liouville measure. Equivalently, 
the latter is the unique semiclassical measure of M. 

The term quantum unique ergodicity refers to the notion of unique ergodicity in 
ergodic theory: a system is uniquely ergodic system if it admits unique invariant 
probability measure. The geodesic flow we are considering admits many invariant 
measures, but the conjecture states that the corresponding quantum system selects 
only one of them. 

3.3.3 Arithmetic quantum unique ergodicity 

The QUE conjecture was motivated by the following result obtained in the cited 
article. Let us recall that any invariant measure can be decomposed into a convex 
combination of ergodic invariant measures: 



/i= / fiedi^ie), (24) 

where {fie, e G Erg} span the set of ergodic probability measure, and i/ is a proba- 
bility measure over this set. The authors specifically considered arithmetic surfaces, 
obtained by quotienting the Poincare disk H by certain congruent co-compact lat- 
tices r. As explained in §2.4.2[ on such a surface it is natural to consider Hecke 
eigenstates, which are joint eigenstates of the Laplacian and the (countably many) 
Hecke operatorq^. It was shown in |90] that for any semiclassical measure fj.sc asso- 
ciated with a sequence of Hecke eigenstates, the ergodic decomposition flMj) of fisc 
does not charge any ergodic component fi^ associated with a periodic geodesic. The 
methods of [90] were refined by Bourgain and Lindenstrauss [27], who showed that 
the measure /i^c of an e-thin tube along any geodesic segment is bounded from above 
by Ce^/^; this bound implies that the Kolmogorov-Sinai entropy^ of almost every 
ergodic component of fisc is bounded from below by 2/9. Finally, using advanced 
ergodic theory methods, Lindenstrauss completed the proof of the QUE conjecture 
in the arithmetic context. 

Theorem 2 [Arithmetic QUE] ^72/ Let M = T\M be an arithmetic surface of con- 
stant negative curvature. Consider an eigenbasis {ipj)j^fis of Hecke eigenstates of Am- 
Then, the only semiclassical measure associated with this sequence is the Liouville 
measure. 



^•^The spectrum of the Laplacian on such a surface is believed to be simple; if this is the case, then an eigenstate 
of A is automatically a Hecke eigenstate. 

'^^The notion of KS entropy will be explained in >jc{.4l 
^^For the precise definition of these surfaces, see |72| . 
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Lindenstrauss and Brooks recently improved this theorem [311 [32]: the QUE result 
holds true, assuming the (ipj) are joint eigenstates of the Laplacian and of a (nonele- 
mentary) single Hecke operator. Their proof uses a new delocalization estimate for 
regular graphs [5U], which yields the same positive entropy results as in [27]. 

This arithmetic QUE result was preceded by a similar statement for the hyper- 
bolic symplectomorphisms on T^ introduced in §2.5.11 These quantum maps Un{So) 
have the nongeneric property to be periodic, so one has an explicit expression for 
their eigenstates. It was shown in [IQ] that for a certain family hyperbolic matrices 
5*0 and certain sequences of prime values of A^, all the eigenstates of Un{So) be- 
come equidistributed in the semiclassical limit, with an explicit bound on the rate of 
equidistribution; they used explicit expressions for the eigenstates, in terms of certain 

2 1 



exponential sums. Some eigenstates, corresponding to the matrix Sdegi — \ o n 

are plotted in Fig. [HI using the Husimi representation on T^. A few years later, 
Kurlberg and Rudnick [65j showed that this equidistribution holds true provided the 
period of the propagator U]\f{So) is large enouglo, which happens to be the case for 
almost all integers N in the limit A^ — ?■ oo. 

A little earlier [61], they had constructed, attached to a given symplectomor- 
phism Sq and any value of N, a finite commutative family of operators {Uj\[{S')S' e 
C{So,N)} including U]y{So), which they called "Hecke operators" by analogy with 
the case of arithmetic surfaces. They then considered specifically the joint ( "Hecke" ) 
eigenbases of this family, and proved QUE in this "arithmetic setting" : 

Theorem 3 [641 Let So G 5*^2 (Z) be a quantizahle symplectic matrix. For each N > 
0, consider a Hecke eigenbasis {4'n,j)j=i,...,n of the quantum map Un{So). Then, for 
any observable f G C°°(T^) and any e > 0, we have 

-l/4+eN 



{^NJ, fNi^Nj) = J fdflL + 0/,e(iV-^/^+^), 

where fi^ is the Liouville (or Lebesgue) measure on T^. 

The eigenstates considered in [iQ] were already particular instances of these Hecke 
eigenstates. 

In view of these positive results, it is tempting to generalize the QUE conjecture 
to other chaotic systems. 

Conjecture 2 [Generalized QUE] Let $^ be an ergodic Hamiltonian flow on some 

energy shell Se- Then, all eigenstates ipn,i of Hn of energies Efij ^ E become equidis- 
tributed when h —> 0. 

Let K be a canonical ergodic map on T^, and (f/Af(fi:))7vgN o?^ associated quantum 
map. Then, all the eigenstates ofUNin) become equidistributed on T^ when N — )■ cxd. 

An intensive numerical study for eigenstates of a Sinai-like billiard was carried on 
by Barnett [15j. It seems to confirm QUE for this system. 

In the next subsection we will exhibit particular systems for which this conjec- 
ture fails. 



^^basically, the period needs to be larger than N^''^^'^ for some e > 0. 



Vol. XIV, 2010 Anatomy of quantum chaotic eigenstates 



19 






ii 


0.5 - 


• J 




>* K 


- 


* 


-0.5 - 



1 -0.5 



1 - 

0.5 - 
- 


1 

4 


- 



^ 
^ 






Figure 8: Husimi functions of various states in Hn (large values=dark regions). Top: 3 (Hecke) 
eigenstates of the quantum cat map Un{Sdegi), for N = 107. Bottom left, center: 2 eigenstates 
of the quantum baker Un{kb) scarred on the period-2 orbit, for iV = 48 and N = 128. Bottom 
right: random state ([5^ for N — 56. (Reprinted from [S2]) 



3.3.4 Counterexamples to QUE for quantum maps 

In this section we will exhibit sequences of eigenstates of certain quantized chaotic 
maps, converging to seniiclassical measures different from fii, thus disproving the 
above conjecture. 

Let us continue our discussion of symplectomorphisms on T^. We recall that 
for any N > 1, the quantized symplectomorphism Un^Sq) is periodic (up to a 
global phase) of period T^ < 3N, so that its eigenvalues are essentially T^r-roots 
of unity. For values of A^ such that Tj^ <t^ N, the spectrum of Un^Sq) is very 
degenerate, in which case imposing the eigenstates to be of Hecke type becomes 
a strong requirement. We already mentioned that, provided the period is not too 
small (namely, T/v S> iV^/^+'^, which is the case for almost all values of A^), then 
QUE holds for any eigenbasis [65]. On the opposite, there exist (sparse) values of 
A^, for which the period can be as small as T/v ~ C logN, so that the eigenspaces 
have huge dimensions ~ C~^A^/logA^. This freedom allowed Faure, De Bievre and 
the author to explicitly construct eigenstates with different localization properties 



Theorem 4 Take Sq G SL2{'L) a (quantizahle) hyperbolic matrix. Then, there exists 

an infinite (sparse) sequence 5 C N such that, for any periodic orbit 7 of ks^, 

one can construct a sequence of eigenstates {iPn)n€S o/f/iv(S'o) associated with the 

semiclassical measure 

1 1 



l^s 



-/i-y + -/^L • 



(25) 



More generally, for any ks^^ -invariant measure fiinv, one can construct sequences of 
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Figure 9: Left: Husimi density of an eigenstate of UNiScat)i half-scarred on a fixed point. Notice 
the hyperbohc structure around the fixed point (Repr. from [4$). Center, right: two eigenstates 
of the Walsh-quantized baker's map (we plot their "Walsh- Husimi density"); the corresponding 
semiclassical measures are purely fractal (Repr. from \2,)- 



eigenstates associated with the semiclassical measure 



This result provided the first counterexample to the generalized QUE conjecture. 
The eigenstates "converging to" |yU^ -|- ^/iL can be called half-scarred. The coeffi- 
cient 1/2 in front of the singular component of /i^ was shown to be optimal |15], a 
phenomenon which was then better understood when considering the entropy of ^sc 
(see g33D. 

Let us briefly explain the construction of eigenstates half-scarred on a fixed 
point 7 = {xo} of ksq- They are obtained by projecting on any eigenspace the 
Gaussian wavepacket (fy^^ (see ( 12T]) ). Each spectral projection can be expressed as 
a linear combination of the evolved states ?7Af(5'o)"v^xo) for ^ ^ [—Tn/'^,Tn/2 — 1]. 
Now, we use the fact that, for A^ in an infinite subsequence 5 C N, the period T^ 
of the operator Un{So) is approximately the double of the Ehrenfest time 



Tp 



logh ^ 
A 



(26) 



(here A is the positive Lyapunov exponent). The above linear combination can 
be split in two components: during the time range n E [—Te/2,Te/2] the states 
Un{So)^^:x_o remain microlocalized at the origin; on the opposite, for times Te/2 < 
\n\ < Te, these states expand along long stretches of stable/unstable manifolds, and 
densely fill the torus. As a result, the sum of these two components is half-localized, 
half-equidistributed. D 

In Fig. [9] (left) we plot the Husimi density associated with one half-localized 
eigenstate of the quantum cat map Uiy{Scat)- 

A nonstandard (Walsh-) quantization of the 2-baker's map was constructed in 
[2], with properties similar to the above quantum cat map. It allows to exhibit semi- 
classical measures ^f^-y+^fJ^L as in the above case, but also purely fractal semiclassical 
measures void of any Liouville component (see Fig. |9]). 

Studying hyperbolic toral symplectomorphisms on T^'^ for d > 2, Kelmer [63] 
identified eigenstates microlocalized on a proper subtorus of dimension > d. He 
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extended his analysis to certain nonlinear perturbations of ksq- Other very explicit 
counterexamples to generalized QUE were constructed in [35], based on interval- 
exchange maps of the interval (such maps are ergodic, but have zero Lyapunov 
exponents) . 

3.3.5 Counterexamples to QUE: bouncing-ball modes of the stadium billiard 

The only counterexample to (generalized) QUE in the case of a chaotic flow concerns 
the stadium billiard, or surfaces with similar dynamical properties. This billiard 
admits a 1-dimensional family of marginally stable periodic orbits, the so-called 
bouncing-ball orbits hitting the horizontal sides of the stadium orthogonally (these 
orbits form a set of Liouville measure zero, so they do not prevent the flow from being 
ergodic). In 1984 Heller [51] had observed that some eigenstates are concentrated 
in the rectangular region (see Fig. [10]). These states were baptized bouncing-ball 
modes, and studied both numerically and theoretically [55] [9]. In particular, the 
relative number of these modes becomes negligible in the limit K — )■ oo, so they 
are still compatible with quantum ergodicity. Yet, no rigorous existence result of 
these modes was known, until Hassell proved that some high-frequency eigenstates 
of (some) stadia indeed have a positive mass on the bouncing ball orbits. To state 
his result, let us parametrize the shape of a stadium billiard by the ratio /3 between 
the length and the height of the rectangle. 

Theorem 5 ^52] For any e > 0, there exists a subset B^ C [1,2] of measure > 
1 — 4e and a number m{e) > such that, for any /3 G B^, the ^-stadium admits a 
semiclassical measure with a weight > m{e) on the bouncing-ball orbits. 

Although the theorem only guarantees that a fraction m(e) of the semiclassical 
measure is localized along the bouncing-ball orbits, the numerical studies suggest 
that some modes are asymptotically fully concentrated on these orbits. Besides, such 
modes are expected to exist for all ratios /3 > 0. 















:;.v*^""'illl)!fl!l!i!ItHii('"»"^::.i 

\' •»,',»«■»•, Mill, •'/*',*,, / 








Figure 10: Two eigenstates of the stadium billiard (/3 = 1). Left (fc = 39.045): the mode has a scar 
along the unstable horizontal orbit. Right (fc — 39.292): the mode is localized in the bouncing-ball 
region. 



3.4 Entropy of the semiclassical measures 

To end this section on the macroscopic properties of eigenfunctions, let us mention 
a recent approach allowing to constrain the semiclassical measures occurring in a 
chaotic system. This approach, initiated by Anantharaman [1], consists in prov- 
ing nontrivial lower bounds for the Kolmogorov-Sinai (KS) entropy of semiclassical 
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measures, using the dispersion property of the Schrodinger equation induced by the 
hyperbohcity of the geodesic flow. The KS entropy is a common tool in classical dy- 
namical systems [60] . To be brief, the entropy H^sif^) of a $*-invariant probability 
measure /i is a nonnegative number which quantifies the information-theoretic com- 
plexity of /^-typical trajectories. As such it does not direcly measure the localization 
of n, but gives some information about it. Here are some relevant properties: 

• the measure /i^ supported on a single periodic orbit 7 has entropy zero. 

• for an Anosov system (flow or diffeomorphism) , the entropy is connected to the 
unstable Jacobian J"(x) (see Fig. [3]) through the Ruelle-Pesin formula: 

V/i invariant, H^silA — l^ogJ^{x.)dfi{x.), (27) 

with equahty iff /i is the Liouville measure. 

• the entropy is an affine function on the set of probability measures: 
HKs{oi^^i + (1 - «)/W2) = aHKsil^i) + (1 - a)HKs{lJ^2)- 

In particular, the invariant measure a/i^ + (1 ~ Ci)f^L of a hyperbolic symplectomor- 
phism S has entropy (1 — a)X, where A is the positive Lyapunov exponent. 

Anantharaman considered the case of geodesic flows on manifolds M of nega- 
tive curvature, see §2.4.21 She proved the following constraint on the semiclassical 
measures of M: 

Theorem 6 /i/ Let (M, g) be a smooth compact riemannian manifold of negative 
sectional curvature. Then there exists c > such that any semiclassical measure figc 
of {M,g) satisfies Hks{i^sc) > c. 

In particular, this result forbids semiclassical measures from being convex combina- 
tions of measures /i-y. supported on periodic geodesies. A more quantitative lower 
bound was obtained in ^ |4j|, related with the instability of the flow. 

Theorem 7 [4lUnder the same assumptions as above, any semiclassical measure 
must satisfy 

HKsif^sc) > j log r dfisc - ^^~^^^^'°- , (28) 

where d = dim M and Amax is the maximal expansion rate of the flow. 

This lower bound was generalized to the case of the Walsh-quantized baker's map 
[2j, and the hyperbolic symplectomorphisms on T^ [251151] . where it takes the form 
Hks{.I^sc) > f- For these maps, the bound is saturated by the half- localized semi- 
classical measures |(/i-y + /i^). 

The bound fl28|) is certainly not optimal in cases of variable curvature. Indeed, 
the right hand side may become negative when the curvature varies too much. A 
more natural lower bound has been obtained by Riviere in two dimensions: 

Theorem 8 /221 WB^Let (M, g) be a compact riemannian surface of nonpositive sec- 
tional curvature. Then any semiclassical measure satisfies 

where A+ is the positive Lyapunov exponent. 
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The same lower bound was obtained by Gutkin for a family of nonsymmetric baker's 
map [48] ; he also showed that the bound is optimal for that system. The lower bound 
(!29|) is also expected to hold for ergodic billiards, like the stadium; it would not 
contradict the existence of semiclassical measures supported on the bouncing ball 
orbits. 

For higher dimensional Anosov systems, one may conjecture PQ E] the lower 
bound 

HKs{f^sc)>^ fhgrdfxsc (30) 

Kelmer's counterexamples [63j show that this bound may be saturated for certain 
Anosov diffeomorphisms on T^*^. In the case of quantized symplectomorphisms on 
T^'^ (nonnecessarily hyperbolic ones). Riviere [8^ recently obtained a lower bound 
of the form 

HxsifJ'sc) > ^ max ( log \(3\ - Amax/2, O) , 

/3eSpec5 

which is sharper than (1251) . but generally weaker than the conjectured bound (150]) 
in the hyperbolic (Anosov) case. 

In view of the Pesin-Ruelle inequality flTTI) . we notice that the QUE conjectured] 
amounts to improving the entropic lower bound fl25]) . or even the conjecture (1301) . 
to HKsif'sc) > J log J" dusc- 

4 Statistical description 

The macroscopic distribution properties described in the previous section give a 
poor description of the eigenstates, compared with our knowledge of eigenmodes of 
integrable systems. At the practical level, one is interested in quantitative properties 
of the eigenmodes at finite values of h. It is also desirable to understand their 
structure at the microscopic scale (the scale of the wavelength i? ~ ^), or at least 
some mesoscopic scale {h '^ R <^ 1). 

The results we will present are of two types. On the one hand, individual eigen- 
functions will be analyzed statistically, e.g. by computing correlation functions or 
value distributions of various representations (position density, Husimi density). On 
the other hand, one can also perform a statistical study of a bunch of eigenfunctions 
(around some large wavevector K), for instance by studying how global indicators of 
localization (e.g. the norms IIV'jUlp) are distributed. We will not attempt to review 
all possible statistical indicators, but only some "popular" ones. 

4.1 Chaotic eigenstates as random states? 

It has been realized quite early that the statistical properties of chaotic eigenstates 
(computed numerically) could be reproduced by considering instead random states. 
The latter are, so far, the best Ansatz we have to describe chaotic eigenstates. Yet, 
one should keep in mind that this Ansatz is of a different type from the WKB Ansatz 
pointwise describing individual eigenstates of integrable systems. By definition, ran- 
dom states only have a chance to capture the statistical properties of the chaotic 
eigenstates, but have no chance to approximate (in L^) any individual eigenf unction. 
This "typicality" of chaotic eigenstates should of course be put in parallel with the 
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typicality of spectral correlations, embodied by the Random matrix conjecture (see 
J.Keating's lecture). 

A major open problem in quantum chaos is to prove this "typicality" of chaotic 
eigenstates. The question seems as open and difficult as the Random Matrix conjec- 
ture. 

4.1.1 Spatial correlations 

Let us now introduce the ensembles of random states. For simplicity we consider 
the Laplacian on a euclidean planar domain Q with chaotic geodesic flow (say, the 
stadium billiard). As in i^, we denote by k'j the eigenvalue of —A corresponding to 
the eigenmode ipj- 

Let us recall some history. Facing the absence of explicit expression for the 
eigenstates, Voros [lUH §7] and Berry [TS] proposed to (brutally) approximate the 
Wigner measures /ij^ of high-frequency eigenstates ipj by the Liouville measure fiL 
on S. We've seen that this approximation is justified by the quantum ergodicity 
theorem, as long as one investigates macroscopic properties of this Wigner function. 
However, Berry [TS] also showed that this approximation provides nontrivial pre- 
dictions for the microscopic correlations of the eigenfunctions ipj. Indeed, a partial 
Fourier transform of the Wigner function leads to the autocorrelation function de- 
scribing the short- distance oscillations of ip. Let us define the correlation function 
by averaging over some distance i?: in 2 dimensions, 

C^A^,r) ''^ ^ / dyr{y-r/2)^{y+r/2) = ^ / dy [ dpe-'^^/'^df,'^{y,p) 

T^-ti J\y-x\<R ^-K J\y~x\<R J 

and take i? to be a mesoscopic scale k~^ <^ R < 1 in order to average over many 
oscillations of ijj. Inserting fii in the place of n^ then provides a simple expression 
for this function in the range < |r| <^ 1: 

c^,R{^,r)^ voi(fi) • ^^^' 

Such a homogeneous and isotropic expression could be expected from our approx- 
imation. Replacing the Wigner distribution by fii suggests that, near each point 
X G fi, the eigenstate ip is an equal mixture of particles of energy k"^ travelling in all 
possible directions. 

4.1.2 A random state Ansatz 

Yet, the Liouville measure fiL is NOT the Wigner distribution of any quantum 
statcj. The next question is thus |1U1] : can one exhibit a family of quantum states, 
the Wigner measures of which resemble ///,? Equivalently, the microscopic correla- 
tions of which behave like (13 ip ? 

To account for these isotropic correlations. Berry proposed the Ansatz of a 
random superposition of isoenergetic plane waves. One form of this Ansatz reads 

9 1/2 ^ 

AandA^) = ( jvVoh^) ) ^( 5Z "j exp(A;% ■ x)) , (32) 



'^^ Characterizing the function on T'R** which are Wigner functions of individual quantum states is a nontrivial 
question. 
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where (nj)j=i,...,Ar are unit vectors randomly distributed on the unit circle, while 
the coefficients (aj)j=i,...,7v are independent identically distributed (i.i.d.) complex 
normal Gaussian random variables. In order to span all possible velocity directions 
(within the uncertainty principle), one should include N p^ k directions n^ The 
normalization ensures that ||'?/'ran(i,fc||L2(Q) ~ 1 with high probability when k ':^ 1. 

Alternatively, one can replace the plane waves in fl32p by circular-symmetric 
waves, namely Bessel functions. The random state then reads 

M 

AandAr,0) = (yo\m-'/^ J2 bmJ\m\ikr)e'^\ (33) 

m=-M 

where the coefficients bm, are i.i.d. complex Gaussian satisfying the symmetry bm = 
b*_^, and M ^ k. Both these random ensembles asymptotically produce the same 
statistical results when /c — )■ oo. 

The random state il^rand,k satisfies the equation (A + k^)ilj = in the interior of 
Q. Furthermore, it satisfies a "local quantum ergodicity" property: for any observable 
f{x,p) supported in the interior of T*r2, the quantum average {iprand,ki fk-^4'rand,k) ~ 
/^l(/) with high probability (more is known about the statistics of these arerages, 
seegXSD. 

The stronger claim is that, in the interior of Q, the local statistical properties of 
4'rand,k, including microscopic ones, should be similar with those of the eigenstates 
ipj with wavevectors kj ^ k. 

The correlation functions of eigenstates of chaotic planar billiards have been 
numerically studied, and compared with the random models, see e.g. [71} [TT] . The 
agreement with f l3T|) is fair for some eigenmodes, but not so good for others; in 
particular the authors observed some anisotropy in the "true" correlation functions, 
which may be related to some form of scarring (see §4.21) . or to the bouncing-ball 
modes of the stadium billiard. 

The value distribution of the random wavefunction fl32|) is Gaussian, and com- 
pares very well with numerical studies of eigenmodes of chaotic billiards [71]. A 
similar analysis has been performed for eigenstates of the Laplacian on a compact 
surface of constant negative curvature [6] . In this geometry the random Ansatz was 
defined in terms of adapted circular hyperbolic waves. The authors checked that 
the coefficients of the individual eigenfunctions, when expanded in these hyperbolic 
waves, are indeed Gaussian distributed. They also checked that the value distribu- 
tion of individual eigenfunctions ipj{x) is Gaussian to a good accuracy, without any 
exceptions. 

4.1.3 On the distribution of quantum averages 

The random state model also predicts the statistical distribution of the quantum 
averages {ipj, fn'ipj), namely the average of the observable / w.r.to the Wigner distri- 
butions, fiYif)- The quantum variance estimate in the proof of Thm. [1] shows that 
the distribution of these averages becomes semiclassically concentrated around the 
classical value HlH)- Using a mixture of semiclassical and random matrix theory 
arguments, Feingold and Peres [16| conjectured that, in the semiclassical limit, the 
quantum averages of eigenstates in a small energy window hkj G [1 — e, 1 + e] should 
be Gaussian distributed, with mean iiiif) and variance related with the classical 
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variance of /. The latter is defined as the integral of the autocorrelation function 
Cfj{t) (see dZD): 

Var,,(/)= [ Cfj{t)dt. 

A more precise semiclassical derivation |13] , using the Gutzwiller trace formula, and 
supported by numerical computations on several chaotic systems, confirmed both the 
Gaussian distribution of the quantum averages, and showed the following connection 
between quantum and classical variances (expressed in semiclassical notations): 

M ^f\ Varc/(/) 

Vars(/) ~ g -— . (34) 

Here g is a symmetry factor {g = 2 in presence of time reversal symmetry, g = I 
otherwise), and Th = 27rhp is the Heisenberg time, where p is the smoothed density 
of states. In the case of the semiclassical Laplacian on a compact surface or a planar 
domain, the mean density p is drawn from Weyl's law (!23|) . so we get 

Vain{f)^ghVaTM)/yom) 

Equivalently, the quantum variance corresponding to wavevectors kj G [K, i^ + 1] is 
predicted to take the value (when K ^ 1): 

M ff\ 9 Varc;(/) 

Successive numerical studies on chaotic euclidean billiards [ini [15] and manifolds or 
billiards of negative curvature [7] globally confirmed this prediction for the quantum 
variance, as well as the Gaussian distribution of the quantum averages at high fre- 
quency. Still, the convergence to this law can be slowed down for billiards admitting 
bouncing-ball eigenmodes, like the stadium billiard [lOj . 

For a generic chaotic system, rigorous semiclassical methods could only prove 
logarithmic upper bounds for the quantum variance [9l] . 

Var,(/)<C/|log;i|. 

Schubert showed that this slow decay can be sharp for certain eigenbases of the 
quantum cat map, in the case of large spectral degeneracies [95] (as we have seen in 
§3.3.4[ such degeneracies are also responsible for the existence exceptionally localized 
eigenstates, so a large variance is not surprising). 

The only systems for which an algebraic decay is proved are of arithmetic nature. 
Luo and Sarnak [75], and then Zhao |11UJ proved that, in the case of the modular 
domain M = S'L2(Z)\E[ (a noncompact, finite volume arithmetic surface for which 
the Laplacian admits many L^ eigenstates), the quantum variance corresponding to 
high-frequency Hecke eigenfunctions is of the form Vari^(/) = — ^: the polynomial 
decay is the same as in fl35l) . but the coefficient B{f) is equal the classical variance 
"decorated" by an extra factor of arithmetic nature. 

More precise results were obtained for quantum symplectomorphisms on the 
2-torus. Kurlberg and Rudnick [67j studied the distribution of quantum averages 
{\/N{^jjnj, fNipNj), j = l,...,iV}, where the (i/^nj) form a Hecke eigenbasis of 
Un{S) (see §3.3.3p . In the semiclassical limit A^ :^ 1, the variance is asymptotically of 
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the form —jp-, with B{f) an "arithmetically decorated" classical variance. They also 
computd the fourth moment of the distribution, which suggests that the latter is not 
Gaussian, but given by a combination of several semicircle laws on [—2, 2] (or Sato- 
Tate distributions). The fifth moment was recently computed by Rosenzweig, and 
shown compatible with this conjecture [88J. The same semicircle law had been shown 
in [nS] to correspond to the asymptotic distribution of the (position) coefficients of 
the Hecke eigenstates, at least for A^ along a subsequence of "split primes" . 

4.1.4 Maxima of eigenfunctions 

Another interesting quantity is the statistics of the maximal values of eigenfunctions, 
that is their L°° norms, or more generally their L^ norms for p G (2, oo] (we always 
assume the eigenfunctions to be L^-normalized) . The maxima belong to the far tail of 
the value distribution, so their behaviour is a priori uncorrelated with the Gaussian 
nature of the latter. 

The random wave model gives the following estimate [S]: for C > large enough. 



^rand,k oo 



\and,k \\2 



< Cy'logk with high probability when A; — )■ oo (36) 



Numerical tests on some euclidean chaotic billiards and a surface of negative curva- 
ture show that this order of magnitude is correct for chaotic eigenstates |8]. Small 
variations were observed between arithmetic/non-arithmetic surfaces of constant 
negative curvature, the sup-norms appearing slightly larger in the arithmetic case, 
but still compatible with f p6|) . For the euclidean billiards, the largest maxima oc- 
cured for states scarred along a periodic orbit (see §4.2p . 

Mathematical results concerning the maxima of eigenstates of generic manifolds 
of negative curvature are scarce. A general upper bound 

m\oo<Ckf-'^^' (37) 

holds for arbitrary compact manifolds [56j , and is saturated in the case of the stan- 
dard spheres. On a manifold of negative curvature, this upper bound can be improved 
by a factor (log kj)~^, taking into account a better bound on the remainder in Weyl's 
law. 

Once again, more precise results have been obtained only in the case of Hecke 
eigenstates on arithmetic manifolds. Iwaniec and Sarnak [57J showed that, for a 
certain family of arithmetic surfaces, the Hecke eigenstates satisfy the bound 



ll^illoo < C,k- 



5/12+e 



5 



and conjecture a bound 

W\oo<aki (38) 

slightly larger than for the random wave model. More recently, Milicevic [77] showed 
that, on certain arithmetic surfaces, a subsequence of Hecke eigenstates satisfies a 
lower bound 

1/2 



-'jiloo 



^--{(i^) '^-«)} 
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thereby violating the random wave model, but still compatible with the conjecture 
(!38l) . The large values are reached on specific CM-points of the surface, of arithmetic 
nature. 

On certain higher dimensional arithmetic manifolds, Rudnick and Sarnak |90] 
had already identified some Hecke eigenstates with larger values, namely 



^ j 11 oo ^ L- kj 



1/2 



A general discussion of this phenomenon appears in the recent work of Milicevic [78] ; 
the author presents a larger family of arithmetic 3-manifolds featuring eigenstates 
with abnormally large values, and conjectures that his list is exhaustive. 



Scat P = + 



Scat P = + 





Figure 11: L°° (left) and L^ (right) norms of the Husimi densities for eigenstates of the quantum 
symplectomorphism UNiSnEGi) (crosses; the dots indicate the states maximally scarred at the 
origin). The data are compared with the values for maximally localized states (Fig. [71 right), la- 
grangian states (Fig. [3 center), random states (Fig. HI bottom right) and the maximally delocalized 
state. (Repr. from [H2). 



4.1.5 Random states on the torus 

In the case of quantized chaotic maps on T^, one can easily setup a model of random 
states mimicking the statistics of eigenstates. The choice is particularly simple when 
the map does not possess any particular symmetry: the ensemble of random states 
in T-Ln is then given by 

1 



^. 



rand,N 



N 



ttiCi. 



(39) 



where {ei)i=i^,„^N is the orthonormal basis f lTSj) of T-Ln, and the (a^) are i.i.d. normal 
complex Gaussian variables. This random ensemble is f/(A^)-invariant, so it takes 
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the same form w.r.to any orthonormal basis of "Hat. For this ensemble, the Husimi 
function Hvrand]v(^) admits, at each point x G T^, an exponential value distribution 
with average 1 + /7v(x), where f^ is exponentially small when N -^ oo. 

This random model, and some variants taking into account symmetries, have 
been used to describe the spatial, but also the phase space distributions of eigenstates 
of quantized chaotic maps. In [S2] , various indicators of localization of the Husimi 
densities (TTSjl ^l have been computed for this random model, and compared with 
numerical results for the eigenstates of the quantized "cat" and baker's maps (see 
Fig. [TT]). The distributions seem compatible with the random state model, except 
for the large deviations of the sup-norms of the Husimi densities, due to eigenstates 
"scarred" at the fixed point at the origin (see Fig. [S]). 

In [^ the Hecke eigenstates of Un{S), expressed in the position basis (e^) as 
in (159]) . were shown to satisfy nontrivial i°° bounda^l: 



'^N,j\ 



sup I Of I <C,N^/^^\ 



For N along a subsequence of "split primes" , the description of the individual Hecke 
eigenstates is much more precise. Their position coefficients are uniformly bounded, 
ll'^AfjIloo < 2, and the value distribution of individual eigenstates {\ipiyj{i/N)\, i = 
0,...,iV — l|}is asymptotically given by the semicircle law on [0, 2], showing that 
these eigenstates are very different from Gaussian random states. In spite of this fact, 
the value distribution of the Husimi function H^^ (x), which, at each x involves 

^ yN position coefficients, appears to be exponential, like in the random state 
model. 

4.2 Scars of periodic orbits 

Around the time the quantum ergodicity theorem was proved, an interesting phe- 
nomenon was observed by Heller in numerical studies on the stadium billiard |54j . 
He noticed that for certain eigenf unctions, the spatial density |'?/'j(a;)P is abnormally 
enhanced along one or several unstable periodic geodesies. He called such an en- 
hancement a scar of the periodic orbit on the eigenstate ipj. See Fig. [10] (left) and 
Fig. [I2] for scars at low and relatively high frequencies. 

This phenomenon was observed to persist at higher and higher frequencies for 
various euclidean billiards |100] [8], but was not detected on manifolds of negative 
curvature [6J. Could scarred states represent counterexamples to quantum unique 
ergodicity? More precise numerics [15] showed that the weights of these enhance- 
ments near the relevant periodic orbit decay in the high-frequency limit, because 
the areas covered by these enhancements decay faster than their intensities. As a 
result, a sequence of scarred states may still become equidistributed in the classical 
limit [SS1[I5]. 

In order to quantitatively characterize the scarring phenomenon, it appeared 
more convenient to switch to phase space representations, in particular the Husimi 
density of the boundary function \l/(g) = dyip{q), which lives on the phase space 
T*dQ of the billiard map [5^ ^^. In this representation, scars were detected as 
enhancements of the Husimi density "Hij, (x) on periodic phase space points x G 
T*dVt of the billiard map (see Fig. H] bottom). 

^®By construction the Husimi densities are L^-normalized. 

-"^^In the chosen normalization, the trivial bound reads lli/'lloo < N^''^ . 
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Figure 12: A high-energy eigenmode of the stadium bilhard (fc w 130). Do you see any scar? 

Similar studies were performed in the case of quantum chaotic maps on the 
torus, hke the baker's map [22] or hyperbohc symplectomorphisms [S2]. The scarred 
state showed in Fig. [8] (left) has the largest value of the Husimi density among 
all eigenstates of Un{Sdegi), but it is nevertheless a Hecke eigenstate, so that its 
Husimi measure should be (macroscopically) close to fii- In this example, the scar- 
ring phenomenon is a microscopic (or mesoscopic) phenomenon, compatible with 
quantum unique ergodicity. This sequence of scarred eigenstates also featured ab- 
normally large values of the Husimi density (see Fig. [Ill left). 

4.2.1 A statistical theory of scars 

Heller first tried to explain the scarring phenomenon using the smoothed local density 
of states 

S^,.,iE) = J2xiE - E,)\{^^,,^P,)\' = {^^,,xiE ~ Hn)^.,) , (40) 

3 

where ip^^ is a Gaussian wavepacket fl?T]) sitting on a point of the periodic orbit, and 
X £ C'^(R) an energy cutoff. The function x is constrained by the fact that this ex- 
pression is estimated through its Fourier transform, namely the time autocorrelation 
function 

t^{^^,,Ulip^,)x{t), (41) 

where x is the ^-Fourier transform of x- Because we can control the evolution of 
(pxo only up to the Ehrenfest time f l26|l . we must take x supported on the interval 
[—Te/2,Te/2], so that x is essentially supported on an interval of width > tj^-jjt- 

Since Uj^ify^f^ comes back to the point xq at each period T, S'^x(-E') has peaks at 
the Bohr-Sommerfeld energies of the orbit, separated by 2'Kh/T from one another. 
However, due to the hyperbolic spreading of the wavepacket, these peaks have widths 
~ Xh/T, where A is the Lyapunov exponent of the orbit. Hence, the peaks can 
only be significant for small enough A, that is weakly unstable orbits. Even if A is 
small, the width Xh/T becomes much larger than the mean level spacing 1/p ~ 
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Ch"^ in the semi classical limit, so that each peak encompasses many eigenvalues. In 
particular, this mechanism can not predict which individual eigenstate will show an 
enhancement at Xq, nor can it predict the value of the enhancements. 

Following Heller's work, Bogomolny j23] and Berry [19] showed that certain 
linear combinations of nearby eigenstates have an "extra intensity" in the spatial 
density (resp. "oscillatory corrections" in the Wigner density) around a certain num- 
ber of closed geodesies. In the semi classical limit, these combinations also involve 
many eigenstates in some energy window. 

A decade later. Heller and Kaplan developed a "nonlinear" theory of scarring, 
which proposes a statistical definition of this phenomenon [5S]. They noticed that, 
given an energy interval / of width ti^ '^ \I\ <^ h, the distribution of the overlaps 
(Kv'x, V'i)!^? Ej G /} depends on the phase point x: if x lies on a (mildly unstable) 
periodic orbit, the distribution is spread between some large values (scarred states) 
and some low values (antiscarred states). On the opposite, if x is a "generic" point, 
the distribution of the overlaps is narrower. 

This remark was made quantitative by defining a stochastic model for the un- 
smoothed local density of states S'x(-E') (that is, taking x i^i POj) to be a delta 
function), as an effective way to take into account the (uncontrolled) long time re- 
currences in the autocorrelation function fj4T|) . According to this model, the overlaps 
|</2x,'0i)P ill ^^ energy window I should behave like random exponential variables, 
of expectation given by the smoothed local density S^^x{E). Hence, if x lies on a short 
periodic orbit, the eigenstates of energies close to the Bohr-Sommerfeld energies 
(where S^^:x_{E) is maximal) statistically have larger overlaps with y^x, while states 
with energies close to the anti-Bohr-Sommerfeld energies statistically have smaller 
overlaps. The concatenation of these exponential random variables with smoothly 
varying expectations produces a non-exponential distribution, with a tail larger than 
the one predicted by Berry's random model. On the opposite, if x is a "generic" 
point, the expectation should not depend on the energy, and the full distribution of 
the |(y9x,V'i)P remains exponential. 

Although not rigorously justified, this statistical definition of scarring gives 
quantitative predictions, and can be viewed as an interesting "dynamical correction" 
of the random state model 



5 Nodal structures 

For a moment let us focus on the case of a Laplacian A^ on a euclidean billiard, 
or a compact riemannian manifold. After having described the "macroscopic skele- 
ton" of the eigenfunctions tpj, namely their semiclassical measures, and the possible 
large values taken by |^j(a;)| near a periodic orbit or elsewhere, we now focus on 
the opposite feature of these eigenfunctions, namely their nodal sets M^^ = {x e 
r2, ijjj{x) = 0}. Since the eigenfunctions ipj can be chosen real, their nodal set is 
a union of hypersurfaces (in 2 dimensions, nodal lines), which sometimes intersect 
each other, or intersect the boundary dQ. This set separates the connected domains 
where ipj{x) has a definite sign, called nodal domains. The nodal set can be viewed 
as a microscopic skeleton of the eigenfunction ipj-. it fully determines the function 
(up to a global factor), and the typical scale separating two nearby hypersurfaces is 
the wavelength k~^ (or hj in the semiclassical formalism). 
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Figure 13: Left: nodal domains of an eigenfunction of the quarter-stadium with k sa 100.5 (do you 
see the boundary of the stadium?). Right: nodal domains of a random state ([33|) with k — 100 
(Repr. from [26]). 



The study of the nodal patterns of eigenfunctions has a long histor y i n mathe- 
matical physics and riemannian geometry. Except for integrable system^lj, we have 
no explicit knowledge of these sets. However, some global properties are known, 
independently on any assumption on the geometry. In 1923 Courant [38] showed 
that, for the Dirichlet Laplacian on a planar domain i7 C M^, the number of nodal 
domains Uj of the j-th eigenstate (counted with multiplicities) satisfies Uj < j, an 
inequality which is an equality in dimension 1. This upper bound was sharpened by 
Pleijel [HI] for high-frequency eigenstates: 

limsup^ < 0.692. 

5.1 Nodal count statistics for chaotic eigenstates 

The specific study of nodal structures of eigenstates for chaotic billiards is more 
recent. Blum, Gnutzmann and Smilansky [22^ seem to be the first authors using 
nodal statistics to differentiate regular from chaotic wavefunctions. They compared 
the nodal count sequence {C,j = -^)j>i for separable vs. chaotic planar domains, and 
observed different statistical behaviours. In the separable case the nodal lines can 
be explicitly computed: they form a "grid" , defined in terms of the two independent 
invariant actions Ji, /2. The distribution of the sequence is peaked near some value 
C,m depending on the geometry. 



Actually, nodal domains are well-understood only for separable systems, a stronger assumption than integrability. 
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Figure 14: Construction of the random-bond percolation model. Left: starting from a square grid 
where tpix) alternatively takes positive (+) and negative (— ) values, a small perturbation 5^{x, y) 
near each crossing creates an avoided crossing. Right: the resulting positive and negative nodal 
domains can be described by drawing bonds (thick/dashed lines) between adjacent sites (reprinted 
from [25]). 



In the chaotic case, the authors found that very few nodal hues intersect each 
other, and conjectured that the sequence should have the same statistics as in the 
case of the random models ( 132|33|) . therefore showing some universality. The nu- 
merical plot of Fig. [13] perfectly illustrates this assertion. In the random model, the 
random variable corresponding to C,j = -^ should be 



N{k)' 



A;> 1 



where z/(A;) is the number of nodal domains in Q for the random state 
N{k) = °^ ' — is the integrated density of states in fi. 



(42) 



and 



Motivated by these observations, Bogomolny and Schmit proposed a heuris- 
tic percolation model to compute the nodal count statistics in random or chaotic 
wavefunctions [25] . Their model starts from a separable eigenfunction, of the form 
ipo{x, y) = cos{kx/\/2) cos{ky/\/2), for which nodal lines form a grid; they perturb 
this function near each intersection, so that each crossing becomes an avoided cross- 
ing (see Fig. [T3j). Although the length of the nodal set is almost unchanged, the 
structure of the nodal domains is drastically modified by this perturbation. Assum- 
ing these local perturbations are uncorrelated, they obtain a representation of the 
nodal domains as clusters of a critical bond percolation model, a well-known model 
in 2-dimensional statistical mechanics. 

The (somewhat amazing) claim made in |25] is that such a perturbation of the 
separable wavefunction has the same nodal count statistics as a random function, 
eventhough the latter is very different from the former in many respects (e.g their 
value distribution). The underlying idea is the following: the uncorrelated local per- 
turbations 6iIj{x, y) instantaneously transform the microscopic square nodal domains 
into mesoscopic (sometimes macroscopic) "fractal" domains, which then do not vary 
too much upon strengthening the perturbation such as to reach a typical random 
state. The high-frequency limit {k — )■ cxd) for the random state (I32|33p corresponds 
to the thermodynamic limit in the percolation model, a limit in which the statisti- 
cal properties of percolation clusters have been much investigated. The nodal count 
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ratio (H3l) counts the number of clusters on a lattice of Ntot = ^N{k) sites. The 
distribution of the number of domains/clusters ^{k) was computed in this limit |25j : 
it is a Gaussian with asymptotic properties 



^'^(^^^ ^-^- 0.0624, ^^^i#l ^^^ 0.0502 . (43) 



>oo 



N{k) ■ ' N{k) 

Nazarov and Sodin [79j have considered random spherical harmonics on the 2-sphere 
(namely, Gaussian random states within each 2n + 1-dimensional eigenspace), and 
proved that when n — )■ oo the distribution of the number of nodal domains of iprand,n 
becomes peaked near the value a n^, for some constant a > 0. Although they compute 
neither the constant a, nor the variance of the distribution, this result indicates that 
the percolation model may indeed correctly predict the nodal count statistics for 
Gaussian random states. 

Remark 2 We now have two levels of modelization. First, the chaotic eigenstates 
are statistically modelled by the random states ^3^S^33\j. Second, the nodal structure 
of random states is modelled by critical percolation. These two conjectures appeal to 
different methods: the second one is a purely statistical problem, while the first one 
belongs to the "chaotic=random" meta- conjecture. 

5.2 Other nodal observables 

According to the Bogomolny-Schmit percolation mode, the statistical distribution 
of the areas of nodal domains has the form V{A) ~ ^~^^^/^^. Of course, this scaling 
can only hold in the mesoscopic range fc~^ ^ ^ ^ 1, since any domain has an area 

> c/e. 

The number Vj of nodal lines of ipj touching the boundary dfl is also an in- 
teresting observable. The random state model (I33p was used in ^22j to predict the 
following properties: 

(Hk)) k^oo^ \dn\_ Var(z>(A;)) fc^oo^ ^ ^^^^ . 



k 277 ' k 

The above expectation value was rigorously proved by Toth and Wigman [98j within 
the following random model: they considered the eigenstates of the Laplacian Aq 
on some planar domain Q, and defined random linear combinations in frequency 
intervals kj G [K, i^ + 1] : 

Aand,K = 5Z "J'^J'' *^^^) 

with the aj i.i.d. normal Gaussians. 

5.3 Macroscopic distribution of the nodal set 

The volume (or length) of the nodal set of eigenf unctions is another interesting 
quantity. A priori, it should be less sensitive to perturbations than the nodal count Uj. 
Several rigorous results have been obtained for this quantity. Donnelly and Fefferman 
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[l2] showed that, for any d- dimensional compact real-analytic manifold, the {d— 1)- 
dimensional volume of the nodal set of any Laplacian eigenstate ijjj satisfies the 
bounds 

C-i kj < Vold_i ^{tPj) <Ckj, 

for some C > depending on the manifold. The statistics of this volume has been 
investigated for various ensembles of random states [El EQl [91] . The average volume 
grows like cm^, with a constant cm > depending on the manifold. Estimates for 
the variance are more difficult to obtain. Berry ^U\ argued that for the 2-dimensional 
random model fl32l) . the variance should be of order log(fc), showing an unusually 
strong concentration property for this random variable. Such a logarithmic variance 
was recently proved by Wigman in the case of random spherical harmonics of the 
2-sphere [TU2] . 

Counting or volume estimates do not provide any information on the spatial 
localization of the nodal set. At the microscopic level, Briining [33] showed that for 
any compact riemannian manifold, the nodal set M{il}j) is "dense" at the scale of 
the wavelength: for some constant C > 0, any ball B{x,C/kj) intersects Nijpj). 

One can also study the "macroscopic distribution" of the zero set, by considering 
the {d — l)-dimensional riemannian measure on //{ipj): 

V/ G C%M), /if ,(/) '^' [ /(x) rf Vol,_i(x) . (45) 

Similarly with the case of the density |?/'j(a;)p or its phase space cousins, the spatial 
distribution of the nodal set can then be described by the weak-* limits of the 

renormalized measures /if = -^ ^^^ in the high-frequency limit. In the case of 

chaotic eigenstates, the following conjecturq^ seems a reasonable "dual" to the 
QUE conjecture: 

Conjecture 3 Let (M, g) be a compact smooth riemannian manifold, with an ergodic 
geodesic flow. Then, for any orthonormal basis {4'j)j>i, the probability measures /i^. 
weak-* converge to the Lebesgue measure on M, in the limit j — t- oo. 

This conjecture is completely open. One slight weakening would be to request that 
the convergence holds for a density 1 subsequence of eigenstates, like in Thm. [T] A 
similar property can be proved in a complex analytic setting (see §5.5p . For M a real 
analytic manifold, eigenfunctions ipj can be analytically continued into holomorphic 
functions ip"^ in some complex neighbourhood M^ of M. For {ipj)j^s a sequence 
of ergodic eigenfunctions (i.e., admitting the Liouville measure as sole semiclassical 
measure), Zelditch has shown |107j that the measures associated with the (complex) 
nodal set of ijjf converge to some limit measure on M*^; however, his result says 
nothing about the distribution of the real zeros (that is, of the set M{ipf) fl M). 

Once more, it is easier to deal with some class of random states than with 
eigenstates. In the case of random spherical harmonics on the sphere, the random 
ensemble is rotation invariant, so for each level n the expectation of the measure 
'"frandn ^^ ^qual to the (normalized) Lebesgue measure. Zelditch |108] generalized 
this result to arbitrary compact riemannian manifolds M, by considering random 



°To my knowledge, this conjecture was first stated by S. Zelditch 
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Figure 15: Stellar representation for 3 eigenstates of the quantum cat map Un{Sdegi)i the Husimi 
densities of which were shown in Fig. [8l top, in a different order. Can you guess the permutation? 



superpositions of eigenstates of the type (jHj): he showed that the expectation of 
l^i'randK conveiges to the Legesgue measure in the hmit K — > oo. 

At the moment the study of nodal sets of random wavefunctions represents hvely 
field of research in probability theory (SUl [TU3j . so the above results are certainly 
not exhaustive. 



5.4 Nodal sets for eigenstates of quantum maps on the torus 

So far we have only considered the nodal set for the eigenfunctions ipj of the Lapla- 
cian, viewed in their spatial representation ip G L^(i7). These eigenfunctions are real 
valued, due to the fact that A is a real operator. At the classical level, this reality 
corresponds to the fact that the classical flow is time reversal invariant. 

In the case of quantized maps on the torus with time reversal symmetry, the 
eigenvectors ipj also have real position components ipj{£/N) = (e^, ipj), i = 0, . . . , N— 
1. Keating, Mezzadri and Monastra JBlj have defined the nodal domains of such 
states as the intervals {ii,ii + 1, ... ,£2 — 1} on which ipjli/N) has a constant 
sign. Using this definition, they computed the exact nodal statistics for a random 
state (!39|) (with real Gaussian coefficients Oj), and numerically showed that these 
statistics are satisfied by eigenstates of a generic quantized Anosov map (namely, a 
perturbation of the symplectomorphism Sdegi)- 

In a further publication [62] , these ideas are extended to chaotic maps on the 4- 
dimensional torus. In this setting, they show that the random state Ansatz directly 
leads to a percolation model on a triangular lattice, with the same thermodynamic 
behaviour as the model of [25]. From this remark, they conjecture that, in the 
appropriate scaling limit, the boudaries of the nodal domains for chaotic eigenstates 
belong to the universality class of SLEg curves, and support this claim by some 
numerics. 



5.5 Husimi nodal sets 

The position representation il){x) is physically natural, but it is not always the most 
appropriate to investigate semiclassical properties: phase space representations of 
the quantum states offer valuable informations, and are more easily compared with 
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invariant sets of the classical flow. We have singled out two such phase space rep- 
resentations: the Wigner function and the Bargmann-Husimi representation. Both 
can be deflned on T*]R'^, but also on the tori T^'^. The Wigner function is real and 
changes sign, so its nodal set is an interesting observable. However, in this section 
we will focus on the nodal sets of the Husimi (or Bargmann) functions. 

The Gaussian wavepackets (12T!) can be appropriately renormalized into states 
(p^ oc v^x depending antiholomorphically on the complex variable z = x ~ ip. Ks a. 
result, the Bargmann function associated with ip G L^(R'^), 

is an entire function of z. In dimension d = 1, the nodal set of Bip{z) is thus a 
discrete set of points in C, which we will denote by Z^. More interestingly, through 
Hadamard's factorization one can essentially recover the Bargmann function Bip 
from its nodal set, and therefore the quantum state if). Assuming ^ Z^, we have 

Bi;{z) = e°^'+'^^+^ n (^ - ^Z-^*) ^^^^^ ' 

leaving only 3 undetermined parameters. Leboeuf and Voros [70] called the set Z^ 
the stellar representation of ip, and proposed to characterize the chaotic eigenstates 
using this representation. This idea is especially appealing in the case of a compact 
phase space like the 2-torus: in that case the Bargmann function Bip of a state 
tpN G T^N is an entire function on C satisfying quasiperiodicity conditions, so that 
its nodal set is Z^-periodic and contains exactly N zeros in each fundamental cell. 
One can then reconstruct the state ■j/'tv from this set of A^ points on T^ (which we 
denote by Zf^): 

BiPi,{z) = e^ W e{z-z,), 



i&Z. 



t2 



where 6{z) is a flxed Jacobi theta function vanishing on Z+iZ, and 7 a normalization 
constant. This stellar representation is exact, minimal {N complex points represent 
ip G 'Hn = C^) and lives in phase space. The conjugation of these three properties 
makes it interesting from a semiclassical point of view. 

In [To] the authors noticed a stark difference between the nodal patterns of 
integrable vs. chaotic eigenstates. In the integrable case, zeros are regularly aligned 
along certain curves, which were identifled as anti-Stokes lines in the complex WKB 
formalism. Namely, the Bargmann function can be approximated by a WKB Ansatz 
similar with ( |20|) with holomorphic phase functions Sj{z), and anti-Stokes lines are 
deflned by the equations Q{Sj{z) — Sk{z)) = in regions where e^^^^^^^^ and e^^'^^^^^'^ 
dominate the other terms. These anti-Stokes lines are sitting at the "antipodes" of 
the lagrangian curve on which the Husimi density is concentrated (for an example 
of such a state, see Fig. [3, center). 

On the opposite, the zeros of chaotic eigenstates appear to be equidistributed 
across the whole torus (see Fig. [T5l) . like the Husimi density itself. This fact was 
checked on other systems, e.g. planar billiards, for which the stellar representations 
of the boundary functions dnipj{x) were investigated in [99j, leading to similar con- 
clusions. This observation was followed by a rigorous statement, which we express by 
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defining (using the same notation as in the previous section) the "stellar measure" 
of a state iI^n £ 'Hn'- 

Theorem 9 l8^] Assume that a sequence of normalized states {ipN € 'Hn)n>i be- 
comes equidistrihuted on T^ in the limit iV — )■ oo (that is, their Husimi measures 
l^i>N ^^Oj^'* converge to the Liouville measure hl). 

Then, the corresponding stellar measures /if also weak-* converge to fii. 

Using the quantum ergodicity theorem (or quantum unique ergodicity when avail- 
able), one deduces that (almost) all sequences of chaotic eigenstates have asymptot- 
ically equidistrihuted Husimi nodal sets. 

This result was proved independently, and in greater generality, by Shiffman 
and Zelditch [96]. It was also extended by Rudnick to the case of holomorphic cusp 
forms on the modular surface [89j. The strategy is to first show that the electrostatic 
potential 

M^^(x) = N-HogH^Jx) = 2N-Hog \B^n{z)\ - 7t\z\^ 

decays (in L^) in the semiclassical limit, and then use the fact that /i^^ = AnAu^^. 
This use of potential theory (specific to the holomorphic setting) explains why such 
a corresponding statement has not been proved yet for the nodal set of real eigen- 
functions (see the discussion at the end of §5.3p . To my knowledge, this result is the 
only rigorous one concerning the stellar representation of chaotic eigenstates. 

In parallel, many studies have been devoted to the statistical properties of stellar 
representation of random states (139|) . which can then be compared with those of 
chaotic eigenstates. Zeros of random entire functions (e.g. random polynomials) have 
a long history in probability theory, see e.g. the recent review [SD], which mentions 
the works of Kac, Littlewood, Offord, Rice. The topic has been revived in the years 
1990 through questions appearing in quantum chaos [BSl EH ED]- 

For a Gaussian ensemble of random states like 039 p . one can explicitly compute 
the n-point correlation functions of the zeros: besides being equidistrihuted, the 
zeros statistically repel each other quadratically on the microscopic scale A^~^/^ (the 
typical distance between nearby zeros), but are uncorrelated at larger distances. 
Such a local repulsion (which imposes a certain rigidity of random nodal sets) holds 
in great generality, showing a form of universality at the microscopic scale [2T1 [SD] . 

The study of |82l| suggested that the localization properties of the Husimi mea- 
sure (e.g. a scar on a periodic orbit) could not be directly visualized in the distribu- 
tion of the few zeros near the scarring orbit, but rather in the collective distribution 
of all zeros. We thus studied in detail the Fourier coefficients of the stellar measures. 

In the case of the random model (l39l) . the variance of the Fourier coefficients could 
be explicitly computed: for fixed k 7^ and A^ ^ 1, the variance is ~ 7r^(^(3)|k|'^/A^'^, 
showing that typical Fourier coefficients are of size ~ N~^^'^. In the case of chaotic 
eigenstates, we conjectured for each fixed k 7^ an absolute upper bound 
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while the equidistribution of Thin. [9] only forces these coefficients to be o(l). We 
also argued that the presence of scars in individual eigenfunctions could be detected 
through the abnomally large values of a few low-k coefficients. 

To finish this section, let us mention a few recent rigorous results concerning 
zeros of random holomorphic functions, after [80] : 

• central limit theorems and large deviation estimates for the "linear statistics" 
//■^(/), with / a fixed test function (including the characteristic function on a 
bounded domain) 

• comparison with other point processes: the Ginibre ensemble obtained by the 
spectrum of random Gaussian complex matrices; ensembles of randomly de- 
formed lattices. Here comes the question of the "best matching" between a 
random set of points with a square lattice. 

• how to use a set of random zeros to partition the plane into cells. 

To my knowledge these properties have not been studied in the framework of Husimi 
functions of chaotic eigenstates. 

To summarize, the stellar representation provides complementary [dual) infor- 
mation to the macroscopic features of the Husimi measures. The very nonlinear 
relation with the wavefunction makes its study difficult, but at the same time inter- 
esting. 
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